The formation rate of semilocal strings 
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We carry out three-dimensional numerical simulations to investigate the formation rate of semilocal 
strings. We find that the back-reaction of the gauge fields on the scalar field evolution is substantial, 
and leads to a significant formation rate in the parameter regime where the semilocal strings are 
classically stable. The formation rate can be as large as a third of that for cosmic strings, depending 
on the model parameters. 
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I. INTRODUCTION 

It is expected that the early Universe underwent a se- 
ries of phase transitions as it cooled down. Typically, 
such transitions lead to the formation of a network of 
defects jjj, which may be of relevance to a number of 
phenomena including structure formation in the Universe 
and baryogenesis || . The best known examples are topo- 
logical defects such as cosmic strings, whose stability is 
guaranteed by the topological structure of the symmetry 
breaking at the phase transition, but it is also possible 
for non-topological defects to form. An example of the 
latter is semilocal strings ||Q ; the semilocal string model 
is simply the Weinberg-Salam model without fermions or 
W bosons, in the limit in which the SU(2) coupling con- 
stant is set to zero. The only parameter in the theory 
is P = irig/m^, the ratio between the scalar and vec- 
tor masses (squared). Its vacuum manifold is the three- 
sphere S 3 , and has no non-contractible loops. Despite 
this Nielsen-Olesen vortices || may form, and are clas- 
sically stable if /3 < 1 Semilocal strings are closely 
related to electroweak strings §, which can be formed 
in the electroweak phase transition, and so understand- 
ing the formation and evolution of these non-topological 
defects is an important task. 

The semilocal model is characterized by the gauge 
fields having insufficient degrees of freedom to be able 
to completely cancel the scalar field gradients, even away 
from the core of any strings which form. While their 
stability depends on the parameter (3, stable semilocal 
strings are stable not only to small perturbations but 
also to semiclassical tunnelling Further, although 
they can, unlike topological strings, come to an end (in 
what is effectively a global monopole) , they will not decay 
by breaking into smaller segments H. On the contrary, 



Hindmarsh jjj) has conjectured that the long-range in- 
teraction between these monopoles should lead to short 
pieces of strings growing into longer ones. 

We have recently shown ]l0| , using a toy model with 
parallel strings, that semilocal strings can be identified 
by studying the pattern of magnetic flux in a simula- 
tion. For topological strings, an estimate of the forma- 
tion rate in systems with planar symmetry is sufficient 
to determine the three-dimensional rate, since topology 
prevents the strings from having an end p!l| . But non- 
topological strings can terminate, with the flux spreading 
out. The closest to a three-dimensional analytic esti- 
mate for semilocal strings is Hindmarsh's calculation || 
of the average magnetic flux through a correlated area 
at j3 = 0, with the conclusion that vortices are rare. For 
the electroweak string, Nagasawa and Yokoyama Q pro- 
posed a technique based on studying the scalar field alone 
and concluded that the initial density would be negligi- 
bly small. However, both approaches neglect the gauge 
fields, which we have found to play a key role [p~0|| . 

In this paper we aim to estimate the density of semilo- 
cal strings at formation. The usual argument for the for- 
mation of (topological) cosmic strings relies on the vac- 
uum manifold having non-contractible loops which can 
force the existence of closed or infinite lines in which the 
Higgs scalar must have zero value, confining the magnetic 
field to these vortex structures where the symmetry has 
been restored. The lack of topology prevents such argu- 
ments being employed for semilocal strings, and numeri- 
cal methods must be employed. This requires two parts 
— a plausible initial configuration where the field config- 
uration captures the essence of a thermal phase transition 
(principally, the existence of a correlation length) and 
secondly dynamical evolution to allow the strings to 'con- 
dense out' and be counted. The numerical simulation of 
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a network of defects is a difficult problem because of the 
large range of scales involved in the dynamics, but it has 
come within the capability of modern supercomputers. It 
is often tacitly assumed that only topological defects are 
sufficiently robust to form in a phase transition through 
the Kibble mechanism ]jj . Given that it is not presently 
known whether semilocal strings form at a comparable 
density to cosmic strings, or with a completely negligible 
density our target is an order-of-magnitude estimate. 

II. THE SIMULATIONS 

We work in flat space-time throughout. The La- 
grangian for the simplest semilocal string model H is 

-iF M ^-^(|0i| 2 + |0 2 | 2 -l) 2 , (I) 

where 4>i and 4*2 are two equally-charged complex scalar 
fields, A^ is a U(l) gauge field and the associated 
gauge field strength. Notice that the gauge coupling and 
the vacuum expectation value of the Higgs have been set 
to one by choosing appropriate units (the inverse vector 
mass, for length, and the symmetry breaking scale, for 
energy). The only remaining parameter in the theory is 
(3 = m 2 /m^, whose value determines the stability of an 
infinitely long, straight, semilocal string with a Nielsen- 
Olesen profile: it is stable for p < 1, neutrally stable 
for P = 1 and unstable for /3 > 1 For /3 = 1 

there is a family of solutions with the same energy and 
different core widths, of which only the semilocal string 
has complete symmetry restoration in the center || . 

We work in temporal gauge Aq — 0. Splitting the 
scalar fields into four real scalars via <f>\ = ip± + iip2, 
4>2 = V J 3 + iipii the- equations of motion are 

~ V 2 Va + P + i>\ + ^ + 1>l ~ 1) ^ (2) 

+,4 2 ^ a + {-l) b (2AV + V.A) i> h = , 

(where b is the complement of a — I h 2, 3 <-> 4 and 
dots are time derivatives) for the scalar fields and 

Ai - V 2 A, + d z V.A + 2U 1 d i f 2 +ip S di i)>A (3) 
+2Ai (i>t + tPI + tP 2 3 + xtf) = , 

for the gauge fields (i = 1,2,3), together with Gauss' 
law, which here is a constraint derived from the gauge 
choice and used to test the stability of the code, 

2 Ui do i>2 + 4> 3 d ^ + diAt - . (4) 

The arrows indicate asymmetric derivatives. 



This system is discretized using a standard staggered 
leapfrog method; however, to reduce its relaxation time 
we also add an ad hoc dissipation term to each equation 
{r/ipi and 77^1 respectively). This is to allow the strings 
to 'condense out' and be identified. In an expanding Uni- 
verse the expansion rate would play such a role, though 
r\ would typically not be constant. We tested a range 
of strengths of dissipation, and checked that it did not 
significantly affect the number densities obtained. Fur- 
ther, our results are always compared to the cosmic string 
case where dissipation would have the same effect. The 
simulations we display later used 77 = 0.5. 

We set up initial conditions for the scalar field by plac- 
ing the field in vacuum with random phases on a subgrid. 
We then iteratively interpolate by bisection onto the full 
grid; after each bisection the field is shifted into the vac- 
uum before the next bisection takes place. 

Having fixed the initial scalar field configuration, we 
need to set the initial gauge field. One possibility is to 
choose the minimum energy configuration on the (fixed) 
scalar field background (this does not mean that the 
gauge fields will cancel all scalar gradients, as there are 
insufficient gauge degrees of freedom). This can only be 
done exactly numerically, as in Ref. jic(] , but in that 
paper it was also shown that an approximate analytic 
minimization ignoring the magnetic flux term performs 
perfectly well: 

Aj(x) = ipidiip 2 - ipzdiipi + ^3<9iV>4 - ^4<3iV>3 • (5) 

All field momenta are set to zero initially. 

These initial conditions are one possible choice out of 
an infinity of possible ways in which one might try to rep- 
resent conditions resembling a thermal phase transition. 
Rather than attempt a highly-accurate description of the 
transition, which is not necessary since we are aiming at 
order-of-magnitude estimates of the formation rate, we 
need only be confident that reasonable changes to the 
initial conditions will not significantly alter the results. 
This we have already tested in two-dimensional simula- 
tions 10 , where we considered a variety of initial condi- 
tions, including ones which may be closer to the sort of 
thermal environment considered in Ref. |f5| . 

III. THE FORMATION RATE 

To minimize any systematic errors in our analysis, we 
always compare our results to the case of cosmic strings, 
which is obtained by simply ignoring one of the (complex) 
scalar fields, setting ^3 = V>4 = 0. This makes the defects 
topological, and the flux tubes formed now map out the 
locations of winding in the scalar field configurations. In 
these field theory simulations of cosmic strings, we can 
follow the early stages of cosmic string network evolution 
by displaying the density of magnetic flux, and we do the 
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FIG. 1. Part of the large simulation, shown at time 60 and 
time 70. Note several joinings of string segments, e.g. two 
separate joinings on the long central string, and the disap- 
pearance of some loops. The different apparent thickness of 
strings is entirely an effect of perspective. 



same with the semilocal string simulations, as was done 
in two dimensions in Ref. JlO| |. 

The starting configurations obtained by our described 
procedure initially yield a complicated mess of flux. How- 
ever, after a few timesteps this resolves itself into loops 
and open segments of string. We observed a clear interac- 
tion between nearby segments which join to form longer 
segments. Fig. 1 shows two time slices from a single large 
/3 = 0.05 simulation, in a 256 3 box, carried out on the 
Cray T3E at NERSC; these are close-ups showing only 
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FIG. 2. This shows the ratio of total string lengths in a 
semilocal and cosmic string simulation, with (3 = 0.05. The 
different lines show different magnetic flux thresholds, from 
bottom to top they are 0.6, 0.55, 0.5, 0.45 and 0.4 times the 
peak flux of a Nielsen-Olesen vortex. 

part of the simulation box.0 We see a collection of short 
string segments and loops; visually this is very differ- 
ent from a cosmic string simulation where strings cannot 
have ends. As time progresses, the short segments either 
disappear or link up to form longer ones. 

We can immediately conclude from these images that 
the formation rate of semilocal strings is not extremely 
close to zero; the fact that flux tubes are observed in 
our simulations implies that the formation rate cannot 
be much smaller than one per correlation volume. 

Nagasawa and Yokoyama jl^] studied the related case 
of electroweak defects and concluded that the initial den- 
sity would be totally negligible. However, this is not 
in contradiction with our results because our semilocal 
strings arise during the evolution due to back-reaction 
on the gauge fields from the scalar field gradients. This 
enables initially short pieces of string to join up to form 
semilocal strings of reasonable size, an effect not included 
in their analysis. Further, the electroweak string resides 
in a part of parameter space where the defects are dy- 
namically unstable, and in this case we find that all the 
flux dissipates soon after the phase transition. 

In order to quantify the formation rate, we compute 
the total length of string in the simulations, always com- 
paring the semilocal string density to that of a cosmic 
string simulation with the same properties (including dis- 
sipation) . We determine the length by setting a magnetic 
flux threshold and computing the fractional volume of 
the box which exceeds it. In Fig. 2, we plot the length 



* Colour images and animations can be found on the WWW at 
cf pa.berkeley . edu/~borrill/def ects/semilocal .html 
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FIG. 3. The ratio of lengths of semilocal and cosmic 
strings. 



of semilocal string relative to the length found in cos- 
mic string simulations, as a function of time and with 
(3 = 0.05. We see that after a transient during which 
the initial tangle of flux sorts itself out, the system set- 
tles down to a reasonable equilibrium. During that ini- 
tial period the ratio of semilocal strings to cosmic ones 
grows, as cosmic ones are there right from the start due 
to topology while the semilocal ones need time to form. 
Even at late times there is a modest upward trend; we 
identify this as being caused by the periodicity of the 
simulation box, which freezes-in any string crossing the 
box, favouring cosmic string annihilation because of their 
higher density. We take the relative densities of semilocal 
and cosmic strings to be that at time 50 in these simula- 
tions. We shall investigate this trend more thoroughly in 
future work using large simulations. There is a modest 
dependence on the choice of flux threshold, and we set it 
at one-half the flux density of a Nielsen-Olesen vortex. 

Fig. 3 shows the ratio of semilocal and cosmic string 
lengths, as a function of the stability parameter (3. These 
results are derived from 700 simulations (50 semilocal 
and 50 cosmic at each of 7 (3 values) carried out in boxes 
of dimension 64 3 using a Sun Ultra II workstation. Al- 
though the initial correlation length of 16 units is a size- 
able fraction of the box size, we are only interested in 
brief evolution to allow the strings to become identifiable 
and so boundary effects are not important. The error 
bars include the statistical spread between simulations, 
and an estimated 25% systematic error from the length 
counting algorithm (see the spread in Fig. 2) and the vis- 
cosity. Those latter uncertainties are the dominant ones. 
Recalling that the formation rate of cosmic strings is es- 
timated to be of order one per correlation volume (0.88 
in Ref . jllj ) , these results are in excellent agreement with 
the two-dimensional results we found in Ref. [|X0}] - They 
show a significant formation rate for low (3, decreasing 
dramatically as (3 — > 1, beyond which there is no evi- 



dence of semilocal string formation. The small amount 
of string seen in some large-/? simulations is an artifact 
of the viscosity; the flux all dissipates if the viscosity is 
turned off, while it persists if (3 < 1. 



IV. CONCLUSIONS 

Our simulations have shown that dynamically stable 
non-topological defects, such as semilocal strings, can 
form with substantial density; this happens even if the 
configuration immediately after the phase transition has 
no symmetry restoration (recall our initial conditions 
place the scalar field in the vacuum everywhere), and 
is due to the back-reaction of the gauge fields on the 
scalars. For (3 — 0.05 we found that the formation rate 
was around 0.3 times that of cosmic strings, while as (3 
is increased towards the stability/instability transition at 
[3 = 1, the density drops to zero. 

Moreover we have observed short segments of string 
joining to form long strings and loops, in agreement 
with Refs. |9|J|] , and collapsing longitudinally, as invoked 
in some baryogenesis scenarios ]l6| . It is important to 
study such phenomena in the related case of electroweak 
strings, where although the monopoles are of finite size 
their cores can become comparable to the separation as 
sin 2 $w - * 1- The details of these dynamics and their 
implications deserve further study. 
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